Contribution of drifting carriers to the Casimir-Lifshitz and Casimir-Polder 
interactions with semiconductor materials 
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We develop a theory for Casimir-Lifshitz and Casimir-Polder interactions with semiconductor 
or insulator surfaces that takes into account charge drift in the bulk material through use of the 
classical Boltzmann equation. We derive frequency-dependent dispersion relations that give the 
usual Lifshitz results for dielectrics as a limiting case and, in the quasi-static limit, coincide with 
those recently computed to account for Debye screening in the thermal Lifshitz force with conducting 
surfaces with small density of carriers. 
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Introduction.- Propagating waves inside semiconduc- 
tors interact with drifting carriers in the bulk material, 
and this is the basis of phenomena such as solid-state 
traveling- wave amplification [l|, Q ■ Typically, an ultra- 
sonic wave or a microwave incident on a semiconductor is 
amplified when the mean drift velocity of carriers exceeds 
the phase velocity of the propagating wave. The theoret- 
ical description of this phenomenon involves Maxwell's 
equations for the electromagnetic field coupled to the 
classical Boltzmann transport equation to describe the 
motion of charged carriers in the bulk semiconductor. 

In principle, the same type of coupling between prop- 
agating waves and drifting carriers is also present for 
quantum vacuum fluctuations of the electromagnetic field 
in the presence of semiconductor boundaries. Hence, 
one should expect that the complete description of the 
Casimir-Lifshitz force between bulk materials and the 
atom-surface Casimir-Polder force Q should take into 
account the possibility of carrier drift when one of the 
surfaces involved is a semiconductor or a conductor with 
small density of carriers. In this limit, the classical Boltz- 
mann equation can be used to determine the dynamic 
equilibrium between a time- and spatially-varying elec- 
tric field and changes in the charge density within the 
material. 

The effect of material properties on quantum vacuum 
forces is encapsulated in the Lifshitz theory through the 
frequency-dependent reflection amplitudes r£j-(i£re) of 
the j-th material boundary. Here p denotes the polariza- 
tion of incoming waves (transverse electric TE or trans- 
verse magnetic TM), k is their transverse momentum, 
and the reflection amplitudes are evaluated at imaginary 
frequencies u = i£ n , where = 2irnkBT / H are the Mat- 
subara frequencies. The Casimir-Lifshitz pressure be- 
tween two plane semi-spaces separated by a vacuum gap 
d is [3 
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where K% = \Jk 2 + £ 2 /c 2 and the prime in the sum over 
n means that a factor 1/2 is to be included for the n = 
term. Assuming that one of the media is dilute, one 
can derive from Eq.([T|) the Casimir-Polder force on an 
atom above a planar surface It has been shown that 
the form of the plate(s) electrical permittivity used to 
compute the reflection amplitudes via Fresnel relations 
vastly alter the magnitude and form of the force in these 
calculations 0]. The effect of carrier drift in the case 
of dynamic fields has not yet been studied in relation to 
Casimir-like forces, and as we show in this Letter, alters 
the form of the field mode equations. 

Recently Pitaevskii [f| has proposed a theory for the 
thermal Lifshitz force between an atom and a conductor 
with a small density of carriers that takes into account 
the penetration of the static component of the fluctuat- 
ing EM field into the conductor. This approach is quasi- 
static, appropriate for the large distance regime of the 
thermal Lifshitz atom-surface interaction, and is essen- 
tially based on the Debye-Hiickel charge screening Q . In 
this static limit, the reflections coefficients are r™(0) = 
and r™(0) = (eoq — k)/(eoq + k). Here eo is the static 
"bare" dielectric constant of the medium (which does 
not take into account the contribution from current car- 
riers), q = \/k 2 + k 2 , and k 2 = 47re 2 no/eo^BT\ where 
— e is the electron charge and no is the (uniform) car- 
rier density [H]. Note that n = l/Ru is the inverse of 
the Debye radius Rb- For good metals the Debye radius 
is very small (on the order of inter-atomic distances), 
while for semiconductors it is much larger (on the order 
of microns or more). This quasi-static calculation for the 
thermal Lifshitz force interpolates between the ideal di- 
electric limit (d <C Rd) and the good conductor limit 
(d 3> Rd)- We further point out that the Debye-Hiickel 
charge screening effect can produce a large correction to 
an electrostatic calibration because a static field can pen- 
etrate a finite distance into the plates, leading to an er- 
ror in the determination of the plate separation On 
the other hand, it can be expected that screening should 
affect dynamic fields as well; however the phenomeno- 
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logical dispersion relation suggested in Q for dynamic 
fields is shown here to be incorrect. In the following 
we will extend Pitaevskii's calculation beyond the quasi- 
static regime, and compute the frequency-dependent TE 
and TM reflection coefficients for semiconductor 

media taking into account carrier drift. 

Field equations .- For an intrinsic semiconductor, the 
densities of carriers and holes are comparable, but the 
dynamics are different. Here we follow the approach 
in [l|, Q, and treat the carriers and holes as dynami- 
cally equivalent, which roughly doubles the charge den- 
sity. This treatment is very accurate in the quasi-static 
limit. Assuming that there is no external applied field 
on the semiconductor, and all fields have a time depen- 
dency of the form e~ lut , Maxwell's equations take the 
form V x E = ifi Lotl, V x H = — ie(w)u>E + J, and 
V - E = — en/e(u>). Here e(w) is the frequency-dependent 
"bare" permittivity of the semiconductor, that does not 
take into account the contribution from current carriers, 
n is the intrinsic carrier density, and is the permeabil- 
ity of vacuum. The carrier current is J = — ertv, where 
v is the mean velocity of carriers. The fact that carriers 
can drift in the semiconductor is modelled by Boltzmann 
transport equation [l], 0| 
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where m is the effective mass of the charge carriers, vt = 
\fh^Tjm is their mean thermal velocity, and r is the 
carrier relaxation time. Linearizing Eq.Q with respect 
to the ac fields which have a factor e _lwt , one can solve 
it for Vn and then use the result in Maxwell's equation 
to derive the fundamental equation for the electric field 
inside the semiconductor 
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Here Co c = lo c /{1 — iur] 
U) c = 47reno/i/e(o>), fi = 



E = [l+i/i e(o;)wi»]V-(V-E) 
(3) 

and D = D/(l — iu>r), where 
erjm is the mobility of carri- 
ers, and D = v\t is the diffusion constant. Note that 
the frequency-dependent ratio lo c /D = 47re 2 no/e(tt>)fcBT 
coincides with k 2 = 1/i? 2 ) in the quasi-static limit. 

TM and TE reflection amplitudes.- Let us assume that 
the semiconductor occupies the semi-space region z < 
and the region z > is vacuum. Eq.Q allows TM 
and TE solutions. For TM modes e y = 0, so that (the 
phase factors e~ tuJt will be omitted from now on) E = 
[e x (z)jt + e z (z)z]e lkx . Substituting this into Eq.([3]) one 
gets two coupled equations, which can be combined into 
two uncoupled fourth-order differential equations for e x 
and e 2 , namely (d 2 — Vt)(^z — Vl) 6 x — and (<9 2 — 



nf,)(<9 2 — r]1)e z = 0, where [l[ 

rj T = k 2 — /ioe(w)w 2 



1 + 1 — 

LO 



k z -i^[i + i— 



D 



LO 



(4) 
(5) 



z — > — oo are e x (z) = 
e nTZ +A' T e VLZ , where we 



The solutions that vanish for 
A T e T ' TZ +A L e^ LZ and e z (z) = A T , 
assume Re t\t and Re r]i to be positive. The amplitudes 
Ai, and At are arbitrary so far, and A' L = —irn^Ai^jk 
and A' T = —ikAx /tjt- The magnetic field inside the 
semiconductor is H = iyATe riTZ e lkx (k 2 — rj^) j hqlotit ■ 

The boundary conditions on the z = interface are 
HE || , E||, Dj^ and continuous. The latter one is au- 
tomatically satisfied for TM modes, while the other ones 
imply E x and H y continuous, and s(lo)E z continuous. 
Imposing these boundary conditions, and using the ex- 
pressions for the fields inside the semiconductor derived 
above, we obtain the reflection amplitude for fields im- 
pinging from the vacuum side, r = (1 — a)/(l + a), with 



2 I 2 

Ul j c 



imaginary frequencies lo 
tude is 



2 I 2 

i£, the TM reflection ampli- 



Expressed along 



„TM 



where y = — 

A TIL 



k' 2 



_ m)v^ r +e/c I -x 

e{it)yjk 2 + £ 2 /c 2 + x 



(6) 



0, so that 
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Plugging this into Eq.Q one gets two 
[<9 2 + /z e(w)w 2 (l + iuj c /Lo + iDk 2 /u)]e x = 



For TE modes 
equations: 

0, and [d 2 — k 2 + /zoe(w)w 2 (l + iuj c /Lo)]e y = ik[l + 
i/j,Qe(oj)LoD]d z e x . The solutions are e x (z) = Ae^ z and 
e y (z) = Be ,1TZ + Ce l3z . Here A and B are constants, (3 2 = 
—i\i^i{uo)Lobr\ 2 L , and C = ikA{3(l+ifi e(Lo)LoD) / \(3 2 — rfc) 
(we assume Re/3 > 0). The magnetic field inside the 
semiconductor is H = (l/ifj,QU))[— [Br\Te r]TZ +C(3e l3z )Sc + 
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Imposing the boundary conditions on the z = inter- 
face, and upon performing the rotation lo — > i£, we get 
the TE reflection amplitude for fields impinging from the 
vacuum side 
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Note that r\\ = k 2 + \e(ig) + 47rcr(^)/^]^ 2 /c 2 (where 
cr(«£) = <Jo/(l + C r ) an d Co = e 2 naT jm arc the ac and dc 
Drudc conductivities, respectively), so Eq.([7]) gives the 
usual Fresnel TE reflection coefficient with account of 
ac Drude conductivity. On the other hand, Eq.([6]) gives 
a modified Fresnel TM reflection coefficient due to the 
presence of Debye-Hiickel screening and charge drift. 

Limiting cases.- Let us study the behavior of the 
frequency-dependent reflection amplitudes we have de- 
rived above for some interesting limiting cases. (a) 
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T = 300K 




FIG. 1: Behavior of the functions used to compute the 

Casimir-Lifshitz free energy and entropy for semiconductor 
materials with account of drifting carriers. The reflections 
coefficients are given by and Q, parameters are for in- 
trinsic Ge (see text), and the distance is set to d = lum. The 
variation with temperature (in the range T = — 300K) of the 
TE function is not perceptible on the scale of the figure. The 
corresponding functions without account of Debye screening 
and carrier drift correspond to the T = OK plots in this figure. 



0, we have rj 2 -, 



k 2 



Quasi-static limit: When £ 

eoC^c/c 2 and 77I sa k 2 + k 2 = q 2 (recall that lo c /D = k 2 
in the quasi-static limit). Here eo = e(0). Therefore, 
X ~ k 2 /q, and from ijHJ) and (JT]) we obtain that the zero- 
frequency limit of the reflection amplitudes is r™(0) = 
and r™(0) = (eoq—k)/(eoq+k), which coincides with the 
prediction of [5j for the reflection coefficients in the quasi- 
static limit. In consequence, we recover the correct ther- 
mal Lifshitz force between an atom and a surface with 
small density of carriers, and the associated crossover 
between good conductors and ideal dielectrics, (b) Ideal 
dielectric limit: In this case the free charge density is 
small (no <C 1), and the discrete charges are quasi-bound, 
making their effective thermal velocity very small. There- 
fore, uj c /D f=a 1/Aq is small (as in the quasi-static limit 
for ideal dielectrics), where U) c ~ 47re 2 no/m^e(i^) and 
D are both small. Consequently n\ « k 2 + e(i£)£ 2 /c 2 , 
?7l ~ £ 2 / v t ~ * 00 ! an< i X ~ Vt- We recover from and 
(0 the usual expressions for the reflection coefficients for 
ideal dielectrics. 

Free energy and entropy. - The Casimir-Lifshitz free en- 
ergy for two parallel planar media is 



p n=0 J y ' 
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where = 2-Kk-oT/h and A is the area of the plates. 
Note that we have allowed for an explicit dependence 
of the reflection coefficients on temperature. In Fig.l 
we plot the behavior of g^ as a function of the imagi- 
nary frequency u) — z£ and transverse momentum k for 



TM and TE polarizations (the corresponding reflection 
amplitudes are obtained from Eqs.© and 0). As an 
example, wc consider the case of two identical media 
made of intrinsic germanium. The permittivity of Ge is 
known to have a weak dependence on temperature, and 
becomes constant as T goes to zero [|| . It can be approx- 
imately fitted with a Sellmeier-type expression e(i£) = 
£oo +^o(eo - £oo)/(£ 2 + w 2 .), with e « 16.2, ps 1.1, 
and u>o ~ 5.0 x 10 15 rad/sec. The intrinsic carrier density 
varies with temperature as no(T) = ^Jn c n v e~ E ' J l 2kliT ^ 
where E g is the energy gap, and n c (n v ) is the effective 
density of states in the conduction (valence) band [§]. 
The relaxation time r has an exponential dependency on 
temperature, and at low temperatures goes linearly in 
T to a non-zero constant Q. Given typical parameters 
of intrinsic semiconductors, uj c and D/£ are both very 
small in the relevant range of frequencies for the Lifshitz 
formula, and then only the n = TM mode is modified 
significantly. The effect of drifting carriers can therefore, 
to very high accuracy, be fully modeled by the Debye- 
Hiickel screening length. 

Our theory for Casimir forces taking into account the 
possibility of carrier drift in intrinsic semiconductor me- 
dia is compatible with Nernst theorem of thermodynam- 
ics, that states that the entropy S = —dE/dT should 
vanish at zero temperature for a system with a nondc- 
generate ground state. Whether the systems we are con- 
sidering here have nondegenerate ground states remains 
open; however, satisfaction of the Nernst theorem pro- 
vides weak evidence for the possible viability of a theo- 
retical model. Following, for example, the technique in 
|icl | . and using the fact that the intrinsic carrier density 
vanishes exponentially as T — > (which in turn implies 
that the derivatives of 0) with respect to £ and to 

exponentially vanish at zero frequency as T — > 0, see 
Fig. 1), it can be shown ll| that our theory predicts 
that the Casimir-Lifshitz entropy verifies S(T = 0) = 0. 
The same is true for the Casimir-Polder entropy when 
the plate is a semiconductor. 

The ratio of the Casimir-Lifshitz free energies for pure 
germanium and pure silicon for various conductivity 
models is shown in Fig. 2, where the increase of the 
energy due to the finite conductivity as compared to the 
bare permittivity is demonstrated for large distances. In 
one case, the theory of drifting carriers (Eas. (|6|7p ) is used 
to model the interaction of the field with the plates; in 
the other, a simple additive term to the bare permittiv- 
ity, 47r<7o/£j is employed in the usual Fresncl reflection 



coefficients [101 ■ For the drifting carriers, when the plate 



separation becomes much larger that the Debye-Hiickel 
screening length, the plates appear as perfect conductors 
for the TM n = mode, while in the case of the additive 
term, the plates appear as perfect conductors for the TM 
n = mode at distances of the order of At = Hc/hsT 
(~ 7/um at T = 300K), independent of the material prop- 
erties. 
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FIG. 2: Ratio of Casimir-Lifshitz free energies at T = 300K 
for intrinsic semiconductor parallel plates for different con- 
ductivity models as follows: Carrier drift with Debye-Huckle 
screening, and a dc conductivity term 4nao/£ added to the 
bare permittivity. Parameters are as follows: For Ge, the 
Debye length is Ro = 0.68/im and the dc conductivity is 
cto = 1/(43 fi cm); For Si, Ru =24^manda = l/(2.3xl0 5 Q 
cm). 



the material. However, the finite temperature correction 
described in Q and its apparent disagreement with ex- 
periment cannot be addressed within the scope of our 
model which does not apply to metals, where the elec- 
tron density is sufficiently large that the electron gas is 
degenerate, so use of the classical Boltzmann equation is 
not warranted [l5j |. 

It is possible to show that the reflection amplitudes 
derived in this work can be interpreted in terms of "non- 
local" dielectric functions (spatial dispersion) [3]. We 
have shown that these effects can be derived from read- 
ily available material properties, and that only the quasi- 
static limit (zero Matsubara frequency TM mode) is rel- 
evant. In the near future we plan to apply these results 
to an ongoing measurement of the Casimir-Lifshitz force 
between pure germanium plates. 

We would like to acknowledge correspondence with 
L.P. Pitaevskii, C. Hcnkel and F. Intravaia, and discus- 
sions with G.L. Klimchitskaya and V.M. Mostepanenko. 



Although this effect has yet to be demonstrated for the 
Casimir-Lifshitz force, the drifting carrier treatment pro- 
vides a way to include a finite conductivity term in the 
Casimir-Polder force which has been measured between 
an atom and a fused silica plate [T3|. For fused silica 
the Debye-Huckel screening length is expected to be ex- 
tremely large due to the low charge concentration and 
the quasi-bound character of charges that are contained 
in a dielectric material. Taking an effective i?D > 1cm is 
certainly not unreasonable, in which case the fused sil- 
ica used in [l| can be treated as a perfect dielectric, as 
has been assumed in the analysis of this experiment. On 
the other hand, including the effect of dc conductivity by 
adding 4irao/£, to the permittivity in the usual Frcsncl 
formulas leads to a increase in the force (by up to nearly 
a factor of two [T3|) at distances of order 10 /^m where 
the experiment was performed, and this disagrees with 
the experimental result. 

Conclusions.- We have shown that treating the finite 
conductivity of a non-degenerate semiconductor (or in- 
sulator) by use of the classical Boltzmann equation in 
conjunction with Maxwell's equation leads to a modi- 
fication of the Casimir-Lifshitz force between such ma- 
terials and provides a way to describe the effects of a 
small conductivity. In particular, for small electric fields 
such that \eE\Rv/kBT <C 1, as expected for Casimir and 
related forces, a standard treatment of adding a term 
Airiao/uj to a "bare" dielectric permittivity is not cor- 
rect for distances less than the Debye-Huckle screening 
length. This is because the current driven by the electric 
field, J = ctE, is counterbalanced by thermal diffusion, 
as modelled through the classical Boltzmann equation. 
Thus, this result represents the dynamic equilibrium be- 
tween a time- varying field and the charge distribution in 
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